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Abstract 

■ There exists a bijection between the set of Latin squares of order n and the set of 

I feasible solutions of the 3-dimensional planar assignment problem {3PAPn). In this 

paper, we prove that, given a Latin square isotopism Q, we can add some linear con- 
straints to the 3PAPn in order to obtain a 1 — 1 correspondence between the new set of 
feasible solutions and the set of Latin squares of order n having 6 in their autotopism 
group. Moreover, we use Grobner bases in order to describe an algorithm that allows 
one to obtain the cardinal of both sets. 



> 

■ Introduction 

' A Latin square of order n is an n x n array with elements chosen from a set of n distinct 

1^ symbols (in this paper, it will be the set [n] = {1,2, ...,n}) such that each symbol occurs 

, precisely once in each row and each column. The set of Latin squares of order n is denoted 

by LS{n). A partial Latin square of order n, is a n x n array with elements chosen from 
a set of n symbols, such that each symbol occurs at most once in each row and in each 
column. The set of partial Latin squares of order n is denoted by PLS{n). 
^ \ The permutation group on [n] is denoted by 5„. Every permutation 5 ^ Sn can be 

^ ' uniquely written as a composition of pairwise disjoint cycles, 6 = Cf 0C2 o...oCn^, where 

for all i G [115], one has Cf = [cf^i ■■■ )^s^ , with cf ^ = mmj{cf j}. The cycle structure 

of 6 is the sequence I5 = (if, I2, 1^), where if is the number of cycles of length i in 6, for all 
i G [n]. Thus, if is the cardinal of the set of fixed points of 6, Fix{5) = {i G [n] \ 6{i) = i}. 
An isotopism of a Latin square L £ LS{n) is a triple = (a,/3,7) £ In = Sn x Sn x Sn- 
Therefore, a, (3 and 7 are permutations of rows, columns and symbols of L, respectively. 
The cycle structure of is the triple (!„, 1^). 

An isotopism which maps L to itself is an autotopism. The possible cycle structures of 
the set of non-trivial autotopisms of Latin squares of order up to 11 were obtained in Q. 
The set of all possible autotopisms of order n is denoted by 21^. The stabilizer subgroup 
of L in %n is its autotopism group 2t(L). Given Q £ 2l„, the set of all Latin squares L 
such that @ G 2l(L) is denoted by LS{@) and the cardinality of LS{Q) is denoted by 
A(0). Specifically, if 0i and ©2 are two autotopisms with the same cycle structure, then 
A(ei) = A(G2). Given 6 G 21^ and P G PLS{n), the number cp = A(G)/|LSp(G)| is 
called the P-coefficient of symmetry ofQ, where LSp{Q) = {L £ LS{@) \ P C L}. 
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Grobner bases were used in to describe an algorithm that allows one to obtain the 
number A(0) in a computational way. This algorithm was implemented in Singular [^| for 
Latin squares of order up to 7 [^. However, after applying it to upper orders, the authors 
have seen that, in order to improve the time of computation, it is convenient to combine 
Grobner bases with some combinatorial tools. In this paper we study, as a possible tool, 
the 1-1 correspondence between LS{n) and the set of feasible solutions of the 3- dimensional 
planar assignment problem (SPAPn) H: 



min 



S.t. 



Y^iei ^ij*: = l.yj e J,k e K. 

^ijk = 1, Vi e /, fc G X. 
T,keK^ijk = 1,V« G I,j G J. 

x,jk G {0,1}, e l,j e J,k e K. 



(1) 



where wijk are real weights and /, J, K are three disjoint n-sets. Thus, any feasible solution 



of the 3PAPn can be considered as a Latin square L 
K = [n] and Xijk 



(lij) G LS{n), by taking I = J 



0, otherwise. 



. The reciprocal is analogous. 



1 Constraints related to an autotopism of a Latin square 

Given a autotopism 6 = (a,/3,7) G 2l„,, let (1)0 be the set of constraints obtained by 
adding to (1) the constraints xijk = Xa{i)p[j)'y{k)^ ^i ^ I,j ^ J,k ^ K. 

Theorem 1.1. There exists a bijection between LS{Q) and the set of feasible solutions 
related to a combinatorial optimization problem having (l)e as the set of constraints. □ 

(l)e is a system of 3n^ + 2n^ equations of degrees 1 and 2, in n'^ variables, which can 
be solved by using Grobner basis. Thus, if we define F{x) = x ■ {x — 1), then the following 
result is verified: 

Corollary 1.2. LS{Q) corresponds to the set of zeros of the ideal I = ( X]ie[n] -^ijA; ~ ^ I 
j,k G [n]) + (EjGM^ijA; - 1 \i,ke [n]) + ( Efce[n] ^ijfc " 1 \ hj ^ N) +{F{xijk) \ hj,k G 
M) + {xijk - Xa{i)p{j)j{k) I i,j,k G [n]) C Q[x] = Q[xiii, □ 

The symmetrical structure of B can be used to reduce the number of variables of the 

[n], if i ^ Fix{a), 1 
Sp, if i G Fix{a). J 

as a set of (hq, — 1^) • n + 1^ • n/3 multi-indices, where Sa = {cfi \ i G [iIq]} and Sp = {cj^ \ 



previous system. To see it, let us consider S'e = < («, j) | i G Sa,j G 



Proposition 1.3 (Falcon and Martin-Morales Q). Let L = (lij) G LS{Q) be such that all 
the triples of the Latin subrectangle Rl = {{i,j,lij) \ {i,j) G S@} of L are known. Then, 
all the triples of L are known. □ 

Let ifQ be a map in the set of variables x = {xm, ...,Xnnn} such that ipe{xijk) = 

{Xijk, if ihj) € Sq, ^ ^ where m = min{/ G [n] \ (a'(i), /3'(j)) G Sq}. 
Xa"^{i)/3"^{j)'y"^{k)j otherwise. 

Theorem 1.4. LS{@) corresponds to the set of zeros of the ideal L' = (X^jg[„] ipe{xijk) — l \ 
j,k G [n]) + (Eje[n] ^ei^ijk) - 1 \ i,k e [n]) + (Efce[n] V^eixijk) - 1 I ij G N) + {xijk \ 
a{i)=i,m=j,-f{k)^k) +{F{xijk)\{i,j)^Se,ke[n]) =((pe(/)) C Q[(^0(x)]. □ 



Now, let P = (pij) € PLS{n) be such that pij = 0, for all and let cp be 
the P-coefficient of symmetry of Q. Thus, we know that A(G) = cp ■ \LSp{@)\ and we can 
calculate \LSp{Q)\ starting from the set of solutions of an algebraic system of polynomial 
equations associated with Q and P. Specifically, we obtain the following algorithm: 

Algorithm 1.5 (Computation of A(0)). 
Input: Q = (a, /3, 7) G 

Ua, the number of cycles of a; 

P € PLS{n) such that pij = 0, for all Sq] 

Cp, the P-coefficient of symmetry of Q. 

Output: A(0), the number of Latin squares having Q as an autotopism; 

^' •= ( J2te[n] feixijk) -l\j,ke [n]) + { Ejg[„] Veixzjk) -l\i,ke [n]) + ( '£ke[n] feix^Jk) - 1 
hJ e W) + {F{x^Jk) I (ij) eSe,k<E [n]); 

r := /' + {xiji - 4^^^ I Pi J / 0, / G [n] ) + {xup,^^ - 5] \ pi,j ^ 0, / g [n] ) + {xi^p^^^ - 5l \ 
Pij 7^ 0, / € [n] ) ; 5 is Kronecker's delta. 

G := Grobner basis of /' with respect to any term ordering; 

Delta := dimQ(Q[(/90(x)]//'); o Delta is the cardinality of V{r) 

RETURN cp ■ Delta; 



2 Number of Latin squares related to Stg and Qtg. 



We have implemented Algorithm 1.5 in a Singular procedure |^ which improves running 
times of @. Moreover, we have obtained the number A(0) corresponding to autotopisms 
of Slg and Slg, as we can see in Table 1. The timing information, measured in seconds, 
has been taken from an Intel Core 2 Duo Processor T5500, 1.66 GHzvj\th Windows Vista 
operating system. 
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Table 1: Number of Latin squares related to Sis and Slg. 
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